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1.0  INTRODUCTION 


This  report  summarizes  the  results  of  two  studies  conducted  to  address  specific  issues 
related  to  the  utility  of  a  specific  high-resolution  multispectral  sensor,  designated  the  F  sensor.  The 
first  study  assessed  the  validity  of  a  sensor  model,  entitled  the  Performance  Prediction  Model 
(PPM),  which  was  provided  by  the  sensor  contractor.  This  was  performed  by  cross-checking 
results  with  an  independent  sensor  model.  In  most  cases,  good  agreement  was  achieved.  The 
second  study  assessed  the  expected  performance  of  a  vital  component  of  the  F  sensor,  a  magnetic 
bearing  system  (MBS)  for  line-of-sight  control  and  stabilization. 

The  report  is  organized  as  follows.  Section  2  provides  the  PPM  assessment  results. 

Section  3  overviews  the  MBS  analysis.  The  appendix  provides  the  details  of  the  MBS  modeling 
effort. 
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2.0  PPM  ASSESSMENT 


2.1  OBJECTIVE 

The  objective  of  this  task  was  to  perform  an  end-to-end  comparison  of  sensor 
performance  predictions  made  using  the  Performance  Prediction  Model  (PPM),  provided  by  the 
government  [1],  and  the  ERIM  Sensor  Performance  (ESP)  model.  Since  an  exhaustive 
comparison  of  the  two  models  was  not  desired,  a  “black  box”  approach  was  utilized  which  (for 
identical  model  inputs)  compared  PPM  and  ESP  model  output  predictions  for  two  key  sensor 
performance  prediction  metrics: 

1)  Ground  Resolved  Distance  (GRD); 

2)  an  image  quality  metric  based  on  the  Image  Interpreter  Rating  Scale  (IIRS)  [2], 

Several  internal  (to  the  “black  box”)  checks  were  also  made  for  key  intermediate  results: 

1)  Signal-to-noise-ratio  (SNR); 

2)  System  Modulation  Transfer  Function  (MTF). 

The  PPM  was  developed  to  model  the  F  sensor.  Approximately  200  sensor  parameters 
were  supplied  which  define  this  sensor  and  are  used  within  the  PPM  to  predict  the  sensor 
performance.  Rather  than  attempt  to  replicate  the  level  of  detail  within  the  PPM  (a  “bottom-up” 
approach),  the  ESP  “top-down”  modeling  approach  started  with  the  most  significant  parameters 
which  affect  sensor  performance  prediction  and  added  detail  only  if  the  performance  predictions 
were  significantly  affected.  The  ESP  model  had  been  used  to  predict  performance  metrics  for  the 
F  sensor  as  part  of  an  prior  modeling  effort,  using  input  data  provided  by  the  F  sensor  contractor. 
The  input  data  currendy  used  by  the  PPM  varied  slightly  from  that  provided  (or  assumed)  for  the 
prior  study.  In  accordance  with  the  “black  box”  approach,  the  following  input  data  was  taken 
from  the  PPM  and  input  to  the  ESP  model: 

Optical  MTF 

Optical  transmission 

Detector  MTF 

Detector  responsiVity 

Comparison  of  modeling  results  obtained  from  the  ESP  and  PPM  models  for  the  visible  band 
will  be  presented  first,  followed  by  a  comparison  of  modeling  results  for  the  medium  wave  IR 
band. 

2.2  VISIBLE  BAND,  160  KM  METEOROLOGICAL  RANGE  COMPARISON 

The  PPM  contains  a  very  high  visibility  sensor  viewing  scenario  (called  a  “Case” )  with 
the  following  pertinent  parameters: 
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160  km  meteorological  range 
14  degree  solar  zenith  angle 
No  atmospheric  turbulence 
20  km  sensor  altitude,  200  m/s  velocity 

Comparisons  between  the  PPM  and  ESP  models  for  the  160  km  case  were  performed  for  SNR, 
MTF,  GRD,  and  IIRS.  The  following  paragraphs  briefly  describe  the  results  of  the 
comparisons. 

2.2.1  Signal-to-Noise  Ratio 

The  ESP  model  was  run  for  the  above  160  km  met  range  case  (using  MODTRAN)  and 
the  signal-to-noise  calculated  at  16  different  surface  (standoff)  ranges.  Here,  signal  (or  contrast 
signal)  is  defined  as  the  difference  signal  obtained  between  a  detector  element  fully  illuminated 
with  a  15%  target  reflectance  and  a  detector  element  fully  illuminated  with  7%  background 
reflectance.  The  PPM  was  run  for  the  same  viewing  scenario  at  0, 100  and  200  km  surface 
range.  Each  surface  range  calculation  required  approximately  5  minutes  of  computing  time, 
prohibiting  running  the  PPM  at  ail  of  the  ranges.  The  resulting  per  pixel  SNR  values  are 
graphed  in  Figure  2-1  and  compare  very  well,  indicating  that  (in  spite  of  the  different 
implementation  details)  consistent  SNR  results  are  obtained.  This  result  indicates  that  basic 
sensor  model  components  (e.g.  aperture  area,  detector  responsivity,  wavelength  region)  and 
scene  radiance  (determined  from  LOWTRAN  or  MODTRAN)  correspond  closely  between  the 
two  models. 

2.2.2  Modulation  Transfer  Function 

See  the  following  23  km  meteorological  range  visibility  for  MTF  comparison  discussion 
and  plots. 

2.2.3  Ground  Resolved  Distance 

The  ground  resolved  distance  (GRD)  calculation  predicts  the  ground  distance  for  any 
given  slant  range,  SR,  corresponding  to  one  cycle  of  the  standard  Air  Force  Tri-bar  target 
predicted  to  be  at  the  visual  detection  limit.  GRD  is  calculated  as  follows: 


GRD 


SR 

^  limit 


(2-1) 


where  fiimit  is  the  limiting  spatial  frequency,  determined  by  the  spatial  frequency  at  which 
MTFifU  •  SNR(SR)  •  =  SNRD  =  3  (2-2) 

J  limit 
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or 


J  limit  J  Nyquist  ' 

whichever  is  less.  SNRD  is  a  detection  threshold  SNR  determined  empirically  to  be  3  for  the 
visual  task  of  just  resolving  the  tii-bar  target.  The  ratio  of  the  limiting  frequency  to  the  Nyquist 
frequency  is  a  SNR-weighting  term  which  accounts  for  the  fact  that  the  eye  spatially  averages 
the  noise  over  the  larger  bars,  effectively  improving  the  signal-to-noise  at  lower  spatial 
frequencies  (more  pixels  per  bar).  The  limiting  spatial  frequency  is  not  permitted  to  exceed  the 
Nyquist  frequency  (i.e.  is  bounded  above)  due  to  the  confounding  effects  of  aliasing  on  the 
visual  system.  The  GRD  metric  combines  both  MTF  and  SNR  characteristics  of  a  sensor  and  is 
an  indicator  of  image  quality. 


Figure  2-1:  SNR  Comparison  for  Visible  Band,  160  km  Met  Range 


For  this  study,  the  ESP  model  incorporated  an  IIRS  prediction  equation  which  includes 
the  above  conventional  prediction  of  the  GRD,  but  with  an  additional  constraint  placed  on  fiimi, 
which  applies  in  cases  of  high  signal-to-noise  when  the  optical  MTF  limits  the  bandwidth  of  the 
sensor  (as  opposed  to  the  detector  MTF  or  Nyquist  condition  limiting  fumit)-  This  constraint 
does  not  permit  fumii  to  increase  to  unduly  large  values  as  SNR  increases  (e.g.  due  to  a  long, 
high  spatial-frequency  tail  in  the  MTF),  but  recognizes  that  the  system  MTF  introduces  a 
limiting  spatial  frequency  for  conditions  of  high  SNR.  The  method  adopted  to  protect  against 
unduly  small  GRD  due  to  a  long  “tail”  in  the  MTF  (with  little  area  under  the  curve) 
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approximates  an  effective  bandwidth  for  the  system  by  the  following  (heuristic)  technique, 
illustrated  in  Figure  2-2. 

1)  The  inflection  point  of  the  system  MTF  curve  was  found. 

2)  A  tangent  to  the  curve  at  that  point  was  constructed. 

3)  The  intersection  of  that  curve  and  the  spatial  frequency  axis,  fsiope,  was  found. 

4)  This  value  was  used  in  the  calculation  of  GRD  if  it  was  less  than  the  Nyquist 


frequency  and  the  limiting  frequency  determined  from  the  product  of  the  MTF  and 
SNR  (described  above). 


Figure  2-2:  Graphical  Depiction  of  the  Determination  of  fumit 

The  Nyquist  limit  also  protects  against  a  long  MTF  tail,  but  the  above  procedure  provides 
protection  for  those  cases  where  the  optics  (or  LOS  jitter)  are  limiting  the  system  MTF  and  a 
long  system  MTF  tail  (with  little  associated  area)  is  present. 

GRD  was  calculated  for  the  same  set  of  surface  ranges  for  the  two  models.  The  results  are 
shown  in  Figure  2-3.  The  PPM  model  indicates  consistently  smaller  GRD  than  the  ESP.  The 
source  of  the  difference  is  due  to  the  modified  GRD  definition  used  in  the  ESP  (i.e.  the  fsiope  term 
was  less  than  Nyquist  or  the  limiting  frequency  determined  from  the  product  of  the  MTF  and 
SNR  in  this  case).  As  can  be  seen  in  Figure  2-2,  a  fairly  large  area  is  present  under  the  MTF 
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curve  beyond  the  fsiope  point  (in  this  particular  case).  A  better  way  to  determine  an  equivalent 
bandwidth  would  be  to  adjust  fsiope  until  the  area  in  the  triangular  region  is  equal  to  the  area 
under  the  system  MTF  curve.  This  limiting  frequency  would  be  insensitive  to  long  MTF  tails 
with  little  associated  area,  but  would  not  be  as  “heavy  handed”  as  the  current  procedure. 


Figure  2-3:  GRD  Comparison  for  Visible  Band,  160  km  Met  Range 


2.2.4  Image  Quality 

For  this  study,  the  ESP  model  calculated  IIRS  using  the  following  equation: 

IIRS  =  11.6-  3.6  *  \og(^^)  (2-4) 

where  GRD  (in  cm)  was  calculated  as  described  above.  The  slope  and  intercept  were 
determined  by  plotting  the  GRD  values  for  the  IIRS  categories  provided  in  the  NATO  Air 
Standard  101  Document  (Reference  1)  as  shown  in  Figure  2-4. 

The  PPM  model  calculates  IIRS  using  an  image  quality  equation  (IQE)  of  the  following 

form: 

IIRS  =  ai  +  a2  *  \og(RER)  .  as  *  H  -  a4  •  log(GSD)  -  (2-5) 
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where  RER  is  the  relative  edge  response,  GSD  is  the  ground  sampled  distance,  and  H  is  an  edge 
overshoot  correction  term.  The  RER  is  a  number  between  0  and  1  equal  to  the  measured  change 
in  contrast  between  two  pixels  centered  on  a  unit  step  in  input  contrast.  It  is  proportional  to  the 
area  under  the  system  MTF  curve.  Better  optical  systems  (with  wider  bandwidth)  will  have  RER 
values  closer  to  unity.  For  a  sensor  with  RER  close  to  unity  and  large  SNR,  image  quality  will 
be  determined  by  the  ground  sampled  distance  (GSD)  which  is  dependent  only  on  detector  pitch 
and  effective  focal  length  of  the  optical  system.  See  the  PPM  documentation  (Reference  2)  for 
precise  definitions  of  the  RER  and  H  terms.  Values  (from  the  PPM)  used  for  the  coefficients 
were  aj  =  11.8;  a2  =  2.71,  3.32;  a3  =  1.48;  a4  =  3.45,  3.32,  a5  =  0.39.  Coefficients  az  and  84  vary 
for  the  along-track  and  cross-track  case. 


IIRS  image  quality  was  predicted  using  the  PPM  IQE  equation  and  the  ESP  equation 
based  on  GRD.  The  results  are  plotted  in  Figure  2-5.  The  agreement  is  within  0.2  IIRS  at  nadir 
and  0.4  IIRS  at  200  km  surface  range,  a  remarkable  agreement,  given  the  different  approach  used 
to  calculate  IIRS  in  the  two  models.  Consistent  with  the  extremely  good  visibility,  the  dashed 
ESP  IIRS  line  indicates  that  image  quality  is  limited  primarily  by  the  system  MTF  (signal-to- 
noise  is  more  than  adequate)  and  Nyquist-limiting  has  not  yet  been  reached  (i.e.  fsiope  was  used  to 
calculate  GRD). 
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Figure  2-5:  IIRS  Comparison  for  Visible  Band,  160  km  Met  Range 

2.3  VISIBLE  BAND,  23  KM  METEOROLOGICAL  RANGE  COMPARISON 

Model  performance  prediction  comparisons  were  made  for  the  visible  band  at  the  more  realistic 
high-visibility  condition  of  23  km  meteorological  range,  74  degree  solar  zenith  angle,  and 
nominal  turbulence  (Rq  =  0.8  m).  Note  that  northern  hemisphere  mid-latitude  solar  zenith  angle 
is  not  less  than  60  degrees  even  at  mid-day  during  winter  months. 

2.3.1  Signal-to-Noise  Ratio 

The  ESP  model  was  used  to  calculate  SNR  for  various  solar  zenith  angles  and  surface 
ranges.  The  results  are  shown  in  Figure  2-6.  SNR  values  were  calculated  using  the  PPM  at  74 
degrees  solar  zenith  angle.  The  results  are  plotted  in  Figure  2-7.  Again,  the  results  compare 
well,  given  slight  differences  in  model  implementation. 

2.3.2  Modulation  Transfer  Function 

System  MTF  contributors  for  a  nadir  viewing  sensor,  as  modeled  using  the  ESP  model,  are 
shown  in  Figure  2-8.  The  system  MTF  is  determined  primarily  by  the  optical  and  detector 
MTF  characteristics  of  the  sensor.  Smear  due  to  the  sweep  motion  of  the  whiskbroom  in  the 
cross-track  direction  is  greatly  reduced  due  to  the  4-phase  clocking  used  with  the  TDI  detector 
(effectively  permitting  charge  transfer  in  V2 IFOV  increments).  The  sensor  will  exhibit  some 
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aliasing,  as  indicated  by  the  Nyquist  frequency.  Line-of-sight  (LOS)  jitter  was  assumed  (by  the 
PPM)  to  be  negligible.  The  ESP  model  did  not  include  jitter  to  permit  direct  comparison  with 
the  PPM  results.  The  MTF  contributors  at  200  km  surface  range  are  shown  in  Figure  2-9.  Since 
the  whiskbroom  angular  scan  rate  is  constant,  there  is  no  increase  in  blue  due  to  a  longer 
integration  time  at  longer  ranges  (as  would  be  observed  with  a  pushbroom  sensor).  The  system 
MTF  calculated  using  each  model  is  shown  in  Figure  2-10  for  comparison.  The  PPM  system 
MTF  is  lower  than  the  ESP  curve  due  to  a  boundary  layer  MTF  degradation,  not  included  in  the 
ESP  modeling.  This  degradation  is  modeled  (in  the  PPM)  as  a  constant  (with  spatial  frequency) 
0.86. 

2.3.3  Ground  Resolved  Distance 

Predicted  GRD  values  (cm)  for  varying  surface  range  and  solar  zenith  angles  were  computed 
using  the  ESP  model.  The  results  are  shown  in  Figure  2-11.  Up  to  80  km  surface  range,  the 
GRD  is  dominated  by  the  system  MTF  and  is  insensitive  to  the  solar  illumination  level.  A 
comparison  of  the  predicted  GRD  for  the  ESP  and  PPM  models  is  shown  in  Figure  2-12.  The 
GRD  agreement  is  better  than  might  be  expected  since  the  MTF  and  SNR  differences  between 
the  models  compensate  in  the  GRD  calculation. 
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Figure  2-9:  MTF  Components  for  Visible  Band,  200  km  Range,  ESP  Model,  23  km  Met  Range 


Figure  2-10:  MTF  Comparison  for  Visible  Band,  Nadir  Viewing,  23  km  Met  Range 
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Figure  2-11:  GRD  for  Selected  Sun  Angles,  ESP  Model,  23  km  Met  Range 


Figure  2-12:  GRD  Comparison  for  Visible  Band,  23  km  Met  Range 


2.3.4  Image  Quality 


Predicted  IIRS  using  the  ESP  model  are  shown  in  Figure  2-13.  At  low  solar  zenith 
angles  the  SNR  does  not  play  a  significant  role  in  determining  IIRS  (MTF  dominates).  At 
longer  surface  ranges  and  higher  zenith  angles,  the  product  of  the  system  MTF  and  SNR  plays  a 
significant  role  (through  the  GRD  calculation).  The  curve  type  indicates  the  dominating  IIRS 
contributor. 
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Figure  2-13:  IIRS  Predictions  for  Visible  Band,  23  km  Met  Range 

IIRS  image  quality  was  predicted  using  the  PPM  IQE  equation  and  the  ESP  equation 
(based  on  GRD).  The  results  are  plotted  in  Figure  2-14.  The  agreement  is  within  0.2  IIRS  at 
nadir  (as  with  the  1 60  km  met  range  case).  At  200  km  surface  range,  the  ESP  model  indicates 
considerable  degradation  due  to  the  decrease  in  SNR  (and  associated  increase  in  GRD).  The 
PPM  results  are  very  similar  to  the  160  km  met  range  case,  since  the  very  small  coefficient  used 
in  the  IQE  equation  in  the  PPM  model  assigns  little  weight  to  SNR.  Since  the  PPM  predictions 
for  GRD  indicate  an  2.25  factor  increase  at  200  km  relative  to  the  160  km  met  range  result,  the 
PPM  coefficient  for  the  SNR  term  is  questionable  as  is  the  IIRS  prediction  at  longer  ranges  (or 
lower  SNR).  For  every  factor  of  2  increase  in  GRD,  IIRS  is  expected  to  decrease  by  one. 
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Figure  2-14:  IIRS  Comparison  for  Visible  Band,  23  km  Met  Range 

2.4  MID-WAVE  INFRARED  (MWIR)  BAND,  23  KM  METEOROLOGICAL  RANGE 
COMPARISON 

A  comparison  of  the  PPM  and  ESP  models  was  also  performed  for  the  MWIR  band. 
Signal-to-noise,  and  MTF  calculations  were  performed  with  both  models  and  the  results 
compared. 

2.4.1  Signal-to-Noise  Ratio 

Signal-to-noise  calculations  for  the  MWIR  band  using  the  two  models  are  shown  in 
Figure  2-15.  As  in  the  visible  band,  the  results  indicate  that  the  conversion  of  scene  radiance  to 
detected  signal  yields  consistent  values  across  the  two  models. 

2.4.2  Modulation  Transfer  Function 

Comparison  of  the  system  MTF  curves  calculated  using  the  two  models  also  indicate 
very  similar  results  as  shown  in  Figure  2-16.  The  MWIR  band  will  also  exhibit  some  aliasing. 


System  MTF 


2.4.3  Minimum  Detectable  Temperature  Difference  (MDT)  and  Ground  Resolved 
Distance 

Since  SNR  and  MTF  results  are  very  similar  between  the  two  models,  performance 
predictors  such  as  MDT  (which  are  based  on  SNR  and  system  MTF)  can  be  expected  to  be 
similar  as  well.  Calculation  of  derivative  performance  measures  such  as  MDT  or  GRD  were 
not  conducted. 

2.5  SENSOR  COVERAGE  CONSIDERATIONS 

The  PPM  calculates  coverage  for  a  given  IIRS  vale  at  the  mid  point  (in  angle)  of  the 
swath  [see  Figure  2-17].  Once  the  far  point  surface  range  has  been  determined,  the  PPM  model 
can  be  ran  a  second  time  at  that  range  to  determine  the  IIRS  value  at  the  far  end  of  the  swath. 
The  ESP  model  has  adopted  coverage  calculations  where  the  desired  minimum  IIRS  is 
specified  at  the  distant  end  of  the  swath.  Given  that  the  two  models  yield  similar  results  for 
IIRS  as  a  function  of  surface  range  (excepting  the  long  ranges  discussed  above),  similar 
coverage  results  can  be  expected.  Since  comparisons  between  the  two  models  would  indicate 
nothing  more  than  agreement  on  the  calculation  of  the  time  available  for  a  cross  scan  (and  thus 
the  scan  width  on  the  ground),  no  coverage  calculations  were  performed. 


Figure  2-17:  Coverage  Swath  Definitions 


2.6  SUMMARY 


Both  the  PPM  and  ESP  models  yield  consistent  results  for  SNR,  system  MTF,  and  GRD 
predictions.  Remarkably  similar  results  are  obtained  for  the  IIRS  image  quality  predictions, 
since  the  methodologies  to  obtain  IIRS  used  by  the  two  modeling  codes  are  different.  It  is 
suspected  that  the  coefficient  for  the  SNR  term  used  in  the  PPM  is  much  too  small,  since  a 
factor  of  two  increase  in  GRD  at  long  surface  range  (and  low  SNR)  is  not  reflected  in  a 
decrease  in  predicted  IIRS.  Since  the  PPM  IQE  is  similar  in  form  to  that  used  by  the  image 
exploitation  community,  this  equation  is  the  preferred  equation  to  use;  the  IQE  coefficients 
must  first  be  reviewed  and  updated  as  necessary. 

Although  very  detailed  in  its  treatment  of  the  F  sensor,  the  PPM  is  tailored  to  the 
whiskbroom  sensor  and  would  require  considerable  effort  to  modify  for  sensors  which  depart 
significantly  from  the  F  sensor.  Considerable  knowledge  and  facility  with  MathCad  and  Excel 
Macro  worksheets  and  their  interactions  must  be  mastered  before  modifications  can  be  made  by 
the  user. 
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3.0  MBS  STABILITY  ANALYSIS  AND  MODELING 


3.1  INTRODUCTION 

This  report  culminates  the  modeling  and  analysis  development  of  the  MBS  magnetic 
controller.  The  work  is  based  on  information  received  from  the  program  office.  This  task  used 
MathCAD  5  Plus  and  Excel  5.0  for  its  development  environment. 

The  MBS  is  represented  in  Figure  3-1  by  nine  components.  Seven  of  the  components  are 
analog  control,  one  is  digital,  and  the  other  is  the  mass  (J)  of  the  device  under  control.  Of  the 
seven  analog,  three  are  electromagnetic  actuators,  a  torque  motor,  and  3  sets  of  feedback 
sensors.  The  actuators  are  theoretically  identical.  They  can  be  represented  by  a  simple  free 
body  diagram.  In  the  real  physical  system  each  actuator  has  two  degrees  of  freedom.  Each 
degree  is  represented  by  an  independent  circuit.  The  actuator  model  described  in  this  paper  is  a 
generic  representation  of  the  theoretical  system  that  is  the  basis  for  the  six  electromagnets 
designed  into  the  MBS.  The  sensors  shown  in  Figure  3-1  serve  as  feedback  for  actuator  control 
and  thus  are  part  of  the  circuit  modeled. 


Figure  3-1:  MBS  Control  Components 


3.2  LIST  OF  REFERENCE  OR  APPLICABLE  DOCUMENTS 

The  following  is  a  list  of  the  documents  received  from  the  program  office.  The  work 
presented  in  this  report  is  extensively  based  on  the  information  contained  in  these  documents. 
With  the  exception  of  the  drawings  listed,  the  documents  fall  into  two  informational  categories 
They  are  either  design  requirements/specification  or  theoretical  modeling  type  of  documents. 
The  mechanical  drawings  are  dimensionless  pictures  of  the  MBS  stmcture.  None  of  the 
documents  represent  real  tested  characteristics  of  the  MBS. 

1.  Critical  Design  Review  (CDR),  Magnetic  Bearing  Subsystem,  dated  23  April  1993. 
Received  19  July  1993. 

2.  Action  items  for  21  July  Mag  Bearing  CDR.  Received  28  July  1993. 

3.  SatCon  company  brochure.  Received  19  July  1993. 

4.  MBS  CDR  Controller  Design,  dated  7  April  1993.  Received  19  July  1993. 

5.  Magnetic  Bearing  Subsystem  (MBS),  Controller  Presentation,  received  21  July  1993. 

6.  Magnetic  Bearing  Subsystem  (MBS)  Program  Management  Review,  dated  18  Febraary 
1993. 

7.  Magnetic  Bearing  Subsystem  (MBS)  Electronics,  Critical  Design  Review,  dated  21  July 
1993. 

8.  Set  of  9  electrical  schematics  "Red"  stamp  dated  20  July  1993,  received  21  July  1993. 
Drawings  on  file  are: 

•  Drawing  #  (Dwg#),  1072904  (1/2  &  2/2) 

.  Dwg#  1072907 

•  Dwg#  1072905 

.  Dwg#  1072909 

•  Dwg#  1072908 

.  Dwg#  1072903 

•  Dwg#  1072901 

.  Dwg#  1072910  (1/2  only). 

9.  Memorandum,  Response  to  S-100  MBS  Questions  of  3  December  1993. 

10.  Response  to  MBS  Questions  of  16  Feb.  94,  Memorandum  #  9770-94-075. 

11.  MBS  Control  Modeling  Ill,  Memorandum  #  9770-94-067. 

12.  Appendix  to  Memorandum  #  9770-94-075. 

13.  Test  and  Evaluation  Master  Plan  (TEMP),  for  the  MBS. 

14.  Documentation  Tree. 

15.  ES264834:  MBS  Engineering  specification. 

16.  Example  of  Modulation  Transfer  Function  MathCAD  Template  (Sanitized). 

17.  Software  TIM,  MBS  Software  Review  4/6/94. 

18.  CIDS  for  MBS  Control  Algorithms  Spec.  #  C172104. 

19.  SDRL  No.  S013,  Soft.  Requirements  Spec,  for  the  MBS-Controller  CSCI  date  2/17/94. 

20.  SDRL  No.  S022  Soft.  Req.  Spec,  for  the  MBS-Controller  CSCI. 

21.  SDRL  No.  S025,  Rev  B,  Soft.  Req.  Spec,  for  the  MBS-Controller  CSCI. 

22.  SLOC  2-17-94,  Language:  ADA,  source  code  listing  of  6  ADA  files. 

23.  ES261468,  Interface  Control  Document  for  the  MBS. 

24.  Roll  Gimbal  with  Laser  Gyro  &  FSM  Control  Loop  Response;  simulation  source  code. 
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25.  Combined  Vertical/Lateral  Actuator  Model  Forward  Actuator;  simulation  source  code. 

26.  Combined  Vertical/Lateral  Actuator  Model  Aft  Actuator;  simulation  source  code. 

27.  Actuator  Model  Equations. 

28.  Results  of  Error  Sensitivity  Simulations;  Excel  spread  sheet  tables. 

29.  Dwg.  #42294A,  Rec.  4/22/94;  3D  line  drawing  with  no  name  or  number  of  the  Stator 
Assembly. 

30.  Dwg.  #1072488,  Rotor  Assembly  Final  Grind. 

31.  Dwg.  #1072490  Front  Plate-Roll  Motor  Assembly;  2  sheets. 

32.  Dwg.  #422948,  Rec.  4/22/94;  3D  line  drawing  showing  stator  with  fairly  good  view  of  the 
caging  mechanism. 

33.  Dwg.  #1072493,  Rotor  Assembly-  Bonding;  2  sheets. 

34  Dwg.  #1072501,  Core  &  Pole  Assy  FWD  Support  Actuator. 

35  Dwg.  #1072507,  Core  &  Coil  Assy  AFT  Support  Actuator;  2  sheets. 

36.  Dwg.  #1072516,  Pole  &  Core  Assembly  Focus  AFT  Pulling. 

37.  Dwg.  #1072517,  Pole  &  Core  Assembly  Focus,  Forward  Pulling. 

3.3  TECHNICAL  DISCUSSION 

Based  on  the  information  obtained  from  the  references,  the  modeling  was  primarily  done 
in  the  frequency  domain  with  some  work  done  in  the  time  domain  .  The  methods  of  analysis 
were  interpretation  of  Nyquist  contour  plots  and  constant  M-Circles.  The  Nyquist  plots  were 
used  to  demonstrate  that  the  controller  could  be  stable  within  a  range  of  coefficients.  The 
constant  M-Circles  were  used  to  show  that  the  system  could  have  sufficient  bandwidth  to  meet 
requirements  as  stated  in  references  #15  and  #23. 

Figure  3-2  is  an  over-simplified  free  body  diagram  of  the  MBS  structure.  Its  intent  is  to 
show  the  cantilever  configuration  and  the  general  direction  of  the  magnetic  forces  generated  by 
the  actuators.  The  focus  actuators  are  not  shown.  Their  design  is  a  scaled  version  of  the 
Forward  and  Aft  actuators. 

The  type  of  modeling  presented  is  Input/  Output  (I/O),  with  Nyquist  plots  and  constant 
M-Circles.  This  method  looks  at  a  system's  transfer  function  and  determines  the  number  of 
poles  and  zeroes.  From  this  a  graphical  plot  (picture),  is  made  that  can  be  used  to  comment  on 
the  relative  stability  of  the  function.  I/O  models  assume  that  all  initial  conditions  are  zero  and 
they  can  he  generated  analytically  or  empirically.  In  this  analysis  all  plots  were  analytical. 

The  magnetic  actuator’s  transfer  functions,  modeled  Geq(s),  were  derived  from  reference 
#’s  4, 7, 11, 15  and  23.  The  forward  vertical/lateral,  aft  verticaLlateral  and  focus  actuators  use 
the  same  form  of  Geq(s).  To  model  a  specific  actuator  a  profile  is  built  containing  the 
coefficients  unique  to  the  physical  device.  For  example,  the  forward  actuator  has  1 1  poles,  aft 
has  19  and  focus  has  5. 

From  the  models  and  respective  charts  presented,  it  can  be  seen  that  the  transfer 
functions  for  the  actuators  are  stable  for  a  wide  selection  of  gain  and  coefficients.  This  does  not 
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take  into  account  any  disturbance  on  these  systems  from  the  base  plane  source.  But  it  does 
show  that  they  can  be  very  stable. 


Pressure! 


Balkhead 

MBS  Stator 


Base  Plane 


Free  Body  Diagram 

Figure  3-2:  Free  Body  Diagram  of  MBS  Structure 

The  MBS  is  put  in  operation  from  the  caged  mode.  In  the  caged  mode,  the  MBS  is 
caged  with  the  nominal  gap  of  3.1mm.  Prior  to  releasing  the  caged  device,  the  actuators  are 
powered  to  the  level  required  to  suspend  the  rotor  at  the  nominal  gap  (set  point).  The  nominal 
gap  is  a  set  point  generated  by  a  command  from  the  MBS  software.  With  this  scheme,  the 
actuators  are  not  required  to  raise  the  rotor  from  the  safety  supports.  This  may  reflect  a  power 
savings.  Not  having  to  do  this  could  lower  the  duty  cycle  requirements  on  the  power  subsystem. 
In  addition  being  able  to  cage  to  the  nominal  set  point  could  widen  the  area  of  controllability. 
This  could  become  a  critical  point.  We  have  a  linear  controller  operating  on  a  nonlinear 
phenomenon  in  a  vibration  environment  described  by  a  set  of  PSDs  (ref.  #15).  As  seen  by  Eq.2 
of  the  MathCAD  templates  (Appendix  A),  the  magnetic  force  applied  to  the  rotor  is  nonlinear 
with  respect  to  both  gap  (Gap),  and  controller  output  current  (I).  Meeting  the  settling,  slew  and 
collision  requirements  could  be  challenging. 

This  modeling  is  useful  in  helping  to  plan,  develop  or  refine  the  MBS  testing  program. 
As  real  test  data  is  published  on  the  current  device,  it  can  be  input  into  the  model  for  correlation. 
As  it  becomes  necessary  to  characterize  the  detailed  effects  of  disturbances  on  the  MBS,  Root 
Locus  and  state  space  representation  models  could  be  investigated.  With  the  addition  of  a 
complete  state  space  model  the  MBS  as  a  system  would  be  fully  characterized.  This  would 


22 


enable  the  modeler  to  investigate  a  whole  range  of  possible  design  configurations  as  the  MBS 
BOM  settles  into  steady  state. 

Figures  3-3  and  3-4  show  the  block  diagrams  of  transfer  functions  Geq(s)  and  Gtot(s), 
respectively.  The  component  transfer  functions  are  defined  in  the  attached  MathCAD  templates. 
Nyquist  contour  plots  were  developed  for  each.  The  parameters  that  form  the  coefficients  are 
the  best  representation  of  the  theoretical  system  at  this  moment.  As  the  latest  design 
documentation  becomes  available,  the  model  could  be  updated. 
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Actuator  Control  System 

Figure  3-3:  Actuator  Control  System  Block  Diagram 

The  Roll  Torque  motor  is  modeled  as  a  second  order  system  with  the  resolver  as  negative 
feedback.  The  motor  has  three  windings  that  are  independently  controlled.  Modeled  are  two 
variations  of  the  motor.  One  with  variable  armature  current  and  the  other  with  variable  field 
current.  The  actual  MBS  torque  motor  is  a  DC  brush  less  motor.  Motors  of  this  type  are 
sometimes  referred  to  as  “pancake”  motors.  At  this  time  not  much  is  known  about  the  real 
characteristics  of  this  motor.  All  that  has  been  mentioned  (verbally),  is  that  it  will  deliver  the 
torque  as  commanded  by  the  SCP.  For  modeling  proposes  this  would  be  a  torque  command  in, 
with  a  torque  out  linear  relationship. 

The  torque  motors  modeled  are  typical  of  the  characteristics  of  a  wide  range  of  off-the- 
shelf  motors.  The  MathCAD  templates  for  typical  torque  motors  is  attached  to  this  report.  The 
MBS  motor  is  a  design  specific  to  the  MBS  using  mature  DC  brush  less  technology.  The 
developers  have  reported  that  their  initial  testing  of  the  motor  has  gone  well.  They  have  tested  it 
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on  both  SL  bread  board  and  MBS  assemblies  with  realistic  loads.  These  test  were  done  with  zero 
base  plane  vibration. 
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Actuator  Control  System 

Figure  3-4:  Actuator  Control  System  Block  Diagram 


3.4  CONCLUSION 

The  MBS  design  has  evolved  significantly  in  the  past  year.  Much  still  needs  to  be  done. 
There  is  currently  a  prototype  built.  It  will  need  to  under  go  laboratory  testing,  debugging  and 
modification  over  the  next  year.  From  a  stability  point  of  view,  more  then  likely  its  main 
challenge  will  be  in  the  design  of  the  flux  and  position  compensators.  These  devices  reside  in 
the  inner  and  outer  control  loops  of  the  hardware  controller.  Without  them  the  MBS  would  not 
be  stable.  As  the  MBS  is  tested  to  verify  its  compliance  with  the  vibration  requirements  it  is 
likely  that  these  devices  could  be  in  the  critical  path.  As  the  MBS’s  response  to  real  loads  is 
determined,  the  compensators  may  need  adjusting  to  mne  into  the  precise  pole  location. 

It  is  recommended  that  the  Air  Force  continue  with  its  modeling  or  simulation  studies. 
Through  these  studies  tremendous  insight  of  the  characteristics  of  the  acquisitioned  technology 
can  be  independently  gained.  Having  this  capability  within  the  Program  Office  would  aid  the 
scheduling  and  budgeting  processes  as  well  as  assist  in  the  application  of  the  technology. 
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APPENDIX 

MBS  MATHCAD  TEMPLATES 


A-l 


Generic  Actuator  Controller  Model 


This  model  can  be  set  up  to  model  either  the  Forward,  Aft  or  Focus  Actuators.  By  varying  the 
coefficients  on  the  numbered  equations,  the  stability  characteristics  of  an  actuator  can  be  studied. 


The  following  variables  are  for  the  definition  of  the  graphics. 

n:=0..  100  Number  of  points  to  be  plotted. 

d;=l  r:=l  Radius  Of  small  arc. 

c  :=0  m  :=0..400  Pole  to  encircled  and  size  of  matrix, 

t  :=.oi-n  Weighted  value  to  scale  plot. 

n 

RR  :=  100  Radius  of  big  arc. 


The  following  equation  are  used  for  plotting  the  contours  of  the  transfer  functions. 
They  define  the  graphics  in  the  complex  plane. 


L(t,x,y.z)  :=z*j  +x-j  -t-Hy- x)-j 
A(t,e,r,d)  :=e-j  +r-exp{j  -(.5- t)-d-n) 


These  equations  draw  the  Nyquist  contours. 

L  is  for  the  vertical  line  and  A  is  for  the  two  arcs. 


P„  :=L(t„,r,RR.c)  P„^ioo  -  A(t„.c.RR.l) 
Pn+.200  :=L(t„.-RR.-r.c)  P„^3oo  := A(t„.c,r.- 1) 


These  equations  calculate  the  matrix  of  points 
to  be  plotted. 


The  following  variables  are  for  the  definition  of  the  electromagnetic  circuit.  These  were 
initially  taken  from  S1 22  source  code  back  in  October  93.  See  reference  document  #12. 
Since  then  they  have  been  varied  to  study  the  sensitivity  of  each  to  the  stability  of  the 


controller. 

O 

II 

r:=100 

inch 0.0254 

F6:=  11.68 

Np:=4 

01d_R;=  36.09 

Ap  ’=  1.0544-inch^ 

Nt:=780 

pO  :=4-re-10r’ 

g0:=3.M0'^ 
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The  following  equations  were  derived  from  reference  document  #'s  4, 11, 12  and  15. 


Lo  ;=nO-Ap-Np- 


(Lm-t-gO) 


Equ.3 


(Np-Nr) 


Equ.4 


Variables  for  the  following  equations. 
lm:=.001  Br  :=1  gap  ;=. 00498  i:=5 


3.^_Br:lm_^_HiO:^  Equ.  5 
(Im+gap)  (Im+gap) 


( Ap-Np-Nt)  ®  compensator  gain  constant. 


B  =0.987  calculated  value 


Kg:=- 


(pO-Ap) 


Equ7,  Feedback  constant. 


2R)pO 

(Ap-Np) 


Equ8. 


BO  =  0.104 


Calculated  value. 
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Ap-B-Np 


Kb  =2.137*10 


KKb;=2-Fb  Equ.10 

KKb  =  23.36  Calculated  value. 


Calculated  value. 


Ma  := 68.  Supported  cantilevered  weight. 


Position  Compensator 


Gp(s) 


.(s+  .6285436744788757886) -(s-f  3.769790402570060698) 


s-(s+ 125.7) 


Flux  Compensator 


Gf(s)  := 


(s-h  1 57.05765407554671968  )-Kv 
s-(s  + 31420.) 


Equ12. 


Actuator  Coil  &  Supported  Mass 


B.1  iMa-s"' 


Equ13. 


G(s)  := 


•Gp(s)-Gf(s) 


Equ14. 
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These  Plots  are  of  P(s)  in  the  complex  plane.  They  are  examples  of  typical  Nyquist  plots.  They  are  drawn  here 
to  verify  that  the  graphics  are  functioning  properly. 


Generic  Nyquist  Contour  Plots  in  the  complex  plane. 


Magnification  of  contour  plot  at  the  origin. 
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The  following  is  the  graphical  definition  of  the  Constant  M-Circles. 


lOn+.OOlmm-O..  100 


10““  M:=2.05 


C(mm)  ;=  '^  --1-  —7^  -exiKl-Tt-j  -.Ol-mm) 

Im^-  1  M^-  1 


Equation  1 6,  is  the  Transfer  function  plotted  for  the  M-Circle. 


n:=0..100 


(0(ii)  :=10-n-t-.001 


Geqq(s) 


G(s)10^ 

(1  -t-Gfafs)-  Gag(s)) 


Equ.16 


Plot  indicating  bandwidth  of  Geq,  'Constant  M-circle'. 


fai(Gcqq(j  -coCn))) 
Im(C(inm))  0 

fa(Gcqq(j  ©(S))) 


Rc(Gcqq(j  •co(n))).Rc(C(mm)) ,Rc(Gcqq(j  -©(S))) 

Fig.  6 
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Math  Model  of  MBS  Roll  Torque  Motors 


The  motors  modeled  here  are  DC  armature  and  field  current  controlled  motors. 
Values  for  the  selected  coefficients  are  typical  of  those  for  off  the  shielf  motors. 
The  load  is  the  MBS  specification. 


The  following  varibles  define  the  graphics  required  to  plot  in  the 
complex  plane. 

R  ;=  100  This  is  the  outer  radius  of  the  contour  plot, 

n  :=0..  100  This  defines  the  numbe  of  points  to  be  generated  during  the  plotting. 

c  := 0  This  defines  the  complex  point  to  circle, 

d  :=  1  This  defines  the  diameter  of  the  small  graphic  circle, 

r  :=  1  This  defines  the  inner  radius  of  the  countour  plot. 

m  ;=  0..  400  This  defines  the  size  of  the  array  to  be  plotted, 

tn  ■-  .01  -n  j^jg  defines  a  weighting  factor  for  the  plot. 


The  following  two  equations  are  for  the  drawing  of  the  vertical  and  arc  lines  of  the 
Nyquist  contour  plot.  They  are  the  same  as  those  shown  on  page1 . 


L(t,x,y,z)  :=z-j  +x-j  +t-(y-x)-j 

^n4-200 


A(t,e,r,d)  :=e-j  +r-exp(j  -(.5- t)'d-7t) 
^n+-ioo 

**04-300 
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Th0  following  coofficisnts  and  ©guations  define  the  Armature  controlled  motor. 


ma  ;=68. 

Mass  of  unit  load  in  Kg.,  from  the  MBS  specification. 

ro  :=.50 

Moment  arm. 

P:=.01 

Viscous  damping  coefficient. 

Tq:=10 

Specified  Torque. 

Kv:=.l 

Back  emf  constant. 

Ra:=40 

Armature  resistance 

Ia:=2 

Armature  current. 

J  :=0.5-ina-ro^ 

Mass  moment. 

La:=4.7-10'^ 

Armature  inductance  in  henry’s 

J=8.5 

Calculated  value. 

II 

Torque  constant. 

K:=5.0 

Arbituary  gain  used  for  investigations. 

=0.001 

J 

Coefficient  calculated  value 

—  =8.511‘10^ 

La 

Coefficient  calculated  value 

^  =125.156 

J-La 

Coefficient  calculated  value 
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Expression  for  an  armature  controlled  DC  motor. 


Nyquist  Contour 


Magified  Origin,  fig.  8 
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I/O  expression  for  the  motor. 


Lai  \J-La  J-La  /  J 


IQ 

JLa 


Equ.  18 


The  gain  K,  was  inserted  to  aid  in  the  analysis  of  this  type  of  motor. 
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Figures  9  and  10  show  that  a  motor  of  this  type  with  the  a  given  set  of  coefficients 
can  be  stable.  This  is  not  a  surprise  conclusion.  Motors  of  this  type  have  been  built 
for  years.  The  reason  for  this  model  is  to  faciliate  the  investigation  of  the  MBS  motor 
when  its  coefficients  are  published. 


Nyquist  Plot  of  Motor  Transfer  function 


Fig.  9 


Fig.  10 


A-13 


TMOTORY.MCD5 


The  following  coefficients  and  equations  define  the  Field  controlled  motor. 


If  :=4.7-10‘^ 

Field  inductance  in  henry's. 

o 

II 

Field  resistance  in  olms. 

If  .-2 

Field  current  in  amps. 

J  =  8.5 

Calculated  mass  moment. 

K:=5.0 

Gain  factor. 

Kt:=^ 

If 

—  =  8.511*10’ 

If 

Calculated  Coefficients 

=0.001 

.J 

—  =  125.156 
JIf 


Field  Controlled  Motor  Transfer  function. 


K- 


Gpf(s)  :=- 


Kt 

JIf 


S'  s-h 


!]fa 


Equ.  20 
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Figures  1 1  and  12  are  the  typical  plot  of  a  stable  Field  controlled  motor. 


Nyquist  Plot  of  Motor  Transfer  function 


Fig.  11 


Fig.  12 
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Model  of  Magnetic  Actuator  Force 


Parameter  values  taken  from  references  #4  &  #5. 


Nt:=780 

Number  of  turns  in  a  couil 

Ho  :=(4.O-7t-M0r’) 

Permeability  of  gap 

Ap  680.0 

Area  of  coil  in 'm  m^ 

Np:=4 

Number  of  poies 

Gap  :=3.1 

Distance  between  rotot  and  stator  in  'mm'. 

i:=0,.1..2.0 

Actuator  driver  current  in  terms  of  Amps. 

lo  := 0.325 

Bias  current  in  Amps. 

Or  if  we  substitute  equation  Eq.21  into  Eq.22. 
F(I,Gap)  :=  K*I>:NE!!2).X  Equ.  23 

\  2  /  Gap2 
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Figure  13  shows  the  non  linear  characterics  of  the  force  function  with  respect  to 
the  drive  current  at  the  nomimal  gap. 


I 

Driver  Cwrent  in  Amps 

Fig.  13 
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